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Consider  five  points  P\,  Pj,  Ps,Pi,  A  independently  and  uniformly  distributed  over  a 
convex  region  K  in  three  dimensional  space.  We  are  interested  in  the  probability  that  the 
five  points  are  the  vertices  of  a  re-entrant  polyhedron,  i.e.  one  of  the  points  falls  within  the 
tetrahedron  formed  by  the  other  four  points.  This  is  a  three-dimensional  generalization  of  the 
Sylvester  problem  which  asks  for  the  probability  that  four  points  independently  and  uniformly 
distributed  over  a  convex  region  in  the  plane  form  a  re-entrant  quadrilateral,  see  Solomon 
(1978). 


Host  insky  (1925)  has  obtained  the  result  for  the  case  when  A  is  a  sphere,  namely  the 
probability  ps  that  five  random  points  in  a  sphere  form  a  re-entrant  polyhedron  is 
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Ps  = 


143' 


Kingman  (1969)  also  obtained  this  result  by  evaluating  the  probability  that  n-f  2  points  falling 
at  random  in  an  n-dimensional  ball  form  a  convex  hull  with  n+2  rather  than  n-f  1  vertices.  It 


is 


For  n=3,  he  shows  the  probability  is 
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However,  explicit  values  of  pjc  for  random  points  in  non-spherical  convex  regions  have 
apparently  not  yet  been  successfully  calculated.  A  simulation  study  is  performed  here  where  the 
probability  p/c  is  approximated  when  IT  is  a  sphere,  a  cube,  a  tetrahedron,  or  an  octahedron. 
Since  we  know  the  exact  answer  for  the  sphere,  an  anchor  point  is  provided  for  simulation  in 
a  sphere. 
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For  easy  programming  purposes,  a  random  point  P  =  (x,  y,  z)  is  generated  in  a  cube  with 
vertices 

(0, 0, 0),  (0, 0, 1),  (0, 1, 0),  (0, 1, 1),  (1, 0, 0),  ( 1, 0, 1),  (1, 1, 0),  (1, 1, 1) 

where  x,  y,  and  z  are  uniformly  distributed  in  [0, 1].  If  K  is  a  sphere,  then  we  are  interested 
only  in  points  belonging  to  the  particular  sphere  with  centre  (0.5, 0.5, 0.5)  and  radius  0.5  .  So 
a  point  generated  inside  the  unit  cube  is  discarded  if 

(x  -  0.5)*  +  {y  -  0.5)*  +  (*  -  0.5)*  >  0.25. 

Since  the  volume  of  this  sphere  is 

Jrrs  a  \  n  0.52 
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and  the  volume  of  the  unit  cube  is  1,  about  48%  of  the  generated  random  points  will  be 
discarded  during  simulation. 

Another  program  employing  spherical  co-ordinates  can  be  produced  as  follows:  generate 
a  point  P  =  (x,  y,  z)  where  x,  y  ,  and  z  are  independently  normally  distributed  with  mean  0 
and  variance  1.  We  then  perform  a  projection  by  letting 

x'  =  _IL_ 

\\P\\ 

J  =  JL 
r  \\P\\ 

✓  »  JL 
in’ll 

where 

in’ll  =  (x*  +  y* +  **)*. 

The  point  F*  =  (x*,/,  z')  is  then  uniformly  distributed  on  the  surface  of  the  sphere  with  centre 
(0,0,0)  and  radius  1  .  Now  generate  u  ~  £7(0, 1)  and  let 

=  vW 
s'-v'S*' 


then  P"  =  (x",  if",  z")  is  uniformly  distributed  in  the  sphere  . 


Of  the  two  programs  described  above  for  the  sphere,  we  employed  the  former.  To  generate 
uniform  random  numbers,  we  used  the  method  of  combining  the  linear  congruential  generator 
and  the  Fibonacci  generator  which  produced  high  resolution,  shuffled  random  numbers  on  [0,1). 
It  is  easy  and  economical  to  generate  uniform  random  numbers  and  so  the  cost  of  discarding 
points  in  the  cube  that  do  not  fall  in  the  convex  region  is  negligible. 

If  the  convex  region  if  is  a  tetrahedron,  then  we  will  consider  only  points  belonging  to 
the  tetrahedron  contained  in  the  unit  cube,  in  particular  that  body  with  vertices 

(0,0,0),(l,0,0),(i,^,0),  and 

This  tetrahedron  has  edge  length  a  =  1  and  volume 

1 

6>/2  6>/2 

so  about  88%  of  the  generated  random  points  will  be  discarded  during  simulation. 


Similarly,  if  if  is  an  octahedron,  then  the  points  considered  are  those  belonging  to  the 
particular  octahedron  with  vertices 


,  1  1  lWn  n  1,  /n  I  lw  1  n  lw  1  1  1.  .  ,  1  1  , 

2  \/2’  l\/2'  *  ’  (°’0,  2  ’  °’  V2’  2),(V2’0,  Vi'  “d  (2>/2’  2>/2’ 


Since  this  octahedron  has  edge  length  a  = 


i 

75 


and  volume 


a*y/2  1 
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about  83%  of  the  generated  random  points  will  be  discarded.  However  since  generating  random 
points  is  cheap,  this  method  is  still  quite  acceptable. 


The  probability  that  the  five  points  Pi,  A>  A»  P\,  and  A  form  a  re-entrant  polyhedron 
is  the  probability  that  one  of  the  points  falls  within  the  tetrahedron  formed  by  the  other  four 
points  .  There  are  five  possible  combinations  of  a  tetrahedron  formed  by  four  points  from  five 
given  points.  We  have  a  re-entrant  polyhedron  if  any  one  of  these  tetrahedra  contains  the  fifth 
point. 


Suppose  we  want  to  check  if  the  point  P%  lies  within  the  tetrahedron  formed  by 
{Pi,  Pi,  Pi,  A),  that  is  we  want  to  check  that  A  and  the  fourth  point  fall  on  the  same  side  of 


the  plane  passing  through  the  other  three  points  .  Consider  the  plane  passing  through 


Pi  = 

P2  -  (xi,yi,Zi) 

Pi-(*i,V s,*a)- 

The  equation  of  this  plane  is 

a*  +  6y  +  e*  =  d  (*) 


where 

a  =  {Vi~  Vi)(*s  -  *1)  ~  (*s  “  *i)(ya  -  |fi) 

b  =  (Zi~  Zt)(Xi  -  Zi)  -  (*S  -  Xt){Zi  -  zj 

c  =  (xa-  xi)(vt  -  y»)  -  (tfs  “  Vi)(xi  ~  xi) 

d  =  ax  i  +  byi  +  tz\. 

By  replacing  the  co-ordinates  of  Pt  and  A  in  the  equation  above,  we  can  obtain  signed  values 
for  Pt  and  Ps  •  If  the  signs  are  the  same  then  we  can  conclude  that  Pt  and  Pj  fall  on  the  same 
side  of  the  plane  formed  by  (Pi, P2, P3)-  The  process  is  repeated  with  the  other  two  planes 
namely  (Pi,P2)Pi)  and  (P3.Ps.P4). 


We  perform  a  simulation  for  each  of  the  cases  when  if  is  a  sphere,  a  cube,  a  tetrahedron, 


or  an  octahedron  . 

Each  run  is  based  on 

100,000  polyhedra,  and  the  results  obtained  are 

displayed  below 

#  of  re-entrant 

proportion  of  re-entrant 

tetrahedra 

tetrahedra 

Sphere 

6274 

0.06274 

Octahedron 

6728 

0.06728 

Cube 

6945 

0.06945 

Tetrahedron 

8604 

0.08604 

The  probability  of  a  re-entrant  polyhedron  in  a  sphere  follows  a  binomial  distribution 


with 


p  =  —  *  0.0629, 
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80  the  approximate  standard  deviation  when  there  are  100,000  trials  is 


/0.0629  x  0.9371  A 

\  - - *  0.00077. 

V  100,000 

Since  we  obtained  p  =  0.06274  which  is  0.00016  away  from  the  true  value  ,  the  simulation  result 
is  quite  plausible.  The  standard  deviation  (or  the  sphere  approximated  from  the  simulations 
is  also  0.00077,  the  same  value  as  that  obtained  using  the  exact  population  proportion  p.  The 
standard  deviations  for  the  tetrahedron,  cube,  and  octahedron  can  be  estimated  by  y  foo.wo 
from  the  simulations  and  are  respectively  0.00089,  0.00080,  and  0.00080.  This  suggests  95% 
confidence  intervals  of  0.0843  •  0.0878  for  the  tetrahedron,  0.0679  •  0.0710  for  the  cube,  and 
0.0657  -  0.0689  for  the  octahedron.  We  also  note  that  the  order  of  the  proportions  agrees  with 
Sylvester’s  conjecture  that  pk  is  minimised  when  if  is  a  sphere  and  maximized  when  K  is  a 
tetrahedron. 

For  the  two-dimensional  problem,  the  probabilities  for  a  re-entrant  quadrilateral  go  from  | 
to  (=0.2955)  as  one  goes  from  the  triangle  to  the  circle.  Thus  for  two  dimensions  we  have 
0.0378  as  the  difference  between  the  maximum  and  minimum  probabilities;  in  three  dimensions 
we  estimate  0.0233  as  the  range.  This  suggests  that  for  higher  dimensions,  a  good  estimate 
for  the  probability  will  not  differ  much  among  convex  regions.  Thus  the  value  obtained  by 
Kingman  for  the  ball  may  serve  well  for  n  >  4  for  all  convex  regions.  These  values  are: 

n  Ps 

4  0.010710192 

5  0.001542677 

6  0.000194894 

7  0.000022118 

8  0.000002293 

9  0.000000220 


Note  that  the  probability  of  a  re-entrant  polyhedron  in  higher  dimensions  is  exceedingly  small. 


-  *  - 
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